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SUMMARY 


Theodorsen’ s propeller theory is used to calculate the performance 
of a dual -rotating propeller having nonideal load dlstrlhutions. The 
application of this theory to either performance or design calculations 
of dual -rotating propellers Is found to Involve a modification to the 
value of the mass coefficient, which is shown to he alfected hy the 
ratio of the spinner radius to the propeller radius. The effect Is 
very Important and, since it has not heen considered in previous work on 
this subject, the required modification and the method of application is 
explained in the present paper. 

By using airfoil data which were obtained from the results of 
special propeller tests, performance calculations are made for an eight- 
blade dual -rotating propeller with flight Mach numbers varying from 0.53 
to 0 . 90 . The calculated results are ccaiQ)ared with experimental results 
for the same propeller. The calculated and experimental power coeffi- 
cients are found to be in good agreement over the entire range investi- 
gated. At a fll^t Mach number of 0.53j "the calculated and experimental 
efficiences are also in close agreement. At flight flach numbers of O .80 
and 0 . 90 , the calculated efficiences are from 3 to 6 percent lower than 
the experimental values. This discrepancy is believed to be due to 
inacciiracies in the drag data. 


INTRODUCTION 


In references 1 to 4 Theodorsen presents the ideal circTolatlon 
f^anctions for single- and dual -rotating propellers as obtained by a 
method of electrical analogy. As a further step the ideal circulation 
functions are integrated to obtain values of the mass coefficient. The 
mass coefficient is used in developing expressions for computing the 
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thrust, energy loss, and the efficiency of ainy propeller^ having the 
ideal distrltution of circulation. The application of these concepts 
to propeller design is described in reference 5* 

The mass coefficient may be interpreted as th.e ratio of the mean 
rearward displacement velocity taken over the entire wake cross section 
to the rearward displacement velocity on the vortex sheets. The mass 
coefficient is used directly in determining the Interference velocities. 
In practical applications there is a large difference in the mass coef- 
ficients of single- and dual -rotating propellers. In view of the fore- 
going interpretation of the mass coefficient, there is also a large 
difference in the interference velocities. Hence, the use of the cir- 
culation functions of single -rotating propellers in the calculation of 
the performance of dual -rotating propellers, as has been done in the 
past, appears inadequate. 

The circulation functions for single-rotating propellers obtained 
in reference 1 are in good agreement with values previously obtained by 
Goldstein and Lock, which have been found to be adequate for use in per- 
formance calculations involving single -rotating propellers at low flight 
Mach numbers. Very little, if any, insist has been gained, however, as 
to the applicability of Theodorsen’s concepts to calculating the perform- 
ance of dual -rotating propellers. 

It is the piorpose of the present paper to test the applicability of 
the dual-circulation functions in performance calculations by direct 
con^iarlson with ea^ierimental results. Experimental results of the same 
propeller for which the calculations are made, an eight -blade 
NACA 3-(3)(05)-05 dual-rotating propeller, are given in reference 6. 

During the course of the present work it was found that the ratio 
of the spinner radius to the propeller radius has an Important effect 
in evaluating the mass coefficient for dual -rotating propellers. In 
propeller design, this effect can greatly influence the propeller dimen- 
sions. The present paper explains tlie role of- the spinner -propeller 
radius ratio in determining the mass coefficient, and gives a method of 
application in the design procedure. 

The fact that the tests of reference 6 cover fli^t speeds well up 
into the transonic range presents an opportunity to test the applica- 
bility to dual rotation of airfoil data recently acquired from the 
result!^ of special propeller tests (reference 7) in the required Mach 
number range. At low speeds these data are, in general, very similar to 
the restilts obtained from standard two-dimensional wind-tunnel tests of 
the same airfoils. As the critical speed is approached and exceeded, 
several differences are found to exist between the two sets of data. 
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In the present paper courpar Isons are made between the calculated 
propeller cliaracteristlcs and the experimental results at flight Mach 
numbers varying from 0.53 to O. 9 O. The airfoil data used for the cal- 
culations at the lower Mach number are In the range where only small 
discrepancies exist between the propeller airfoil data of reference 7 
and standard two-dimensional data. The differences In the airfoil data 
are much more pronounced at the hl^er flight Mach numbers. 


SYMBOLS 


B 

b 





D 

d 

h 

J 

K(x) 

I 


M 




number of propeller blades 
chord of propeller-blade element 
section drag coefficient 
section lift coefficient 
section design lift coefficient' 
Ideal power coefficient 
power coefficient 
torque coefficient 
thrust coefficient 
diameter of propeller 
drag of propeller section 
blade-section maximum thickness 
advance ratio ('V/nl)) 
circulation function 
lift of propeller section 
Mach number of adVEince 


(p/pn3l)5) 

(Q/pn2D5) 


section resultant Mach number 




NACA EM L5IAI7 


propeller rotational speed, rps 

power 

torque 

tip radius 

radius to any blade element 
thrust of propeller 
velocity of advance 

resultant velocity at blade section 
velocity vector + (itnDx)^^ 

rearward displacement velocity of helical vortex surface 

ratio of displacement velocity to forward velocity (w/V) 

radial location of blade element (r/R) 

spinner radius ratio 

angle of attack, degrees 

blade angle, degrees 

tan‘^ c^/cj 


mass coefficient 




K(x)x dx 


effective mass coefficient |2 



1.0 


K(x)x dx 


- (t* 

axial energy-loss factor 

axial energy-loss factor based on ft* 

propeller efficiency (TV/P or JCt/Cp) 
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p mass density of air 

a solidity (Bb/2«r) 

ac^ propeller -element load coefficient 

0 aerodynamic helix angle 

geometric helix angle 
Subscripts: 

F front 

R rear 

0.7R at O.T radius 

INFLUENCE OF SPINHER-HiOEELLER RADIUS RATIO ON THE 
DESIGN OF DUAL -ROTATING EROPET.TERS 



Although the application of Theodor sen's propeller theory to the 
design of dual -rotating propellers Is described In reference the 
influence of the spinner -propeller radius ratio on the mass coefficient 
is not considered. The necessity for Including this ratio as a design 
parameter arises from the large values of the circulation functions 
occurring over the inner radii for dual rotation, as contrasted to the 
very small values ■which occur In single rotation. 

The mass coefficient 



Includes the entire propeller -disc area from x = 0 to x = 1.0, 

■vdiereas. In actual practice, some of the Inne r radii are covered by a 
spinner. There Is no circulation within the spinner area and the spinner 
surface itself may cause Interference effects on the remaining flow 
field. For cases vdiere this Interference effect Is small, as in the 
present case, the mass coefficient can be evaluated as 


pl.O 

k' = 2 J 3K(x)dx 


CD 
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The Interference velocity is obtained by solving the equation 

Pc = 2 kw(1 + w) ^1 + -^ w] (2) 

where Pq is known for the design conditions (reference 5)- Tbe proper 
determination of w, however, requires a modification of equation (2) to 

Pc = 2k '^?(1 + w) (l + ^ w) (3) 

Comparison of equations (2) and (3) shows that, for a given Pc, equa- 
tion (3) yields a larger value of w than does equation (2) because k' 
is less than k . The proper evaluation of the quantity w is very 
Important because this quantity is used directly in calculating the 
element loading coefficients crcj and the aerodynamic helix angles. A 
dual -rotating propeller which is designed neglecting the spinner propel- 
ler radius ratio will tend to operate with heavier than the desired 
loadings at the design conditions with a consequent loss in efficiency. 

The quantity €‘/k' in equation (3) can' be obtained in the manner 
shown in reference 5 substituting k’ for k. For many practical 
purposes, it sho\ild be sufficiently ancurata to use either c/k or k* 
for c’/k' because the Influence of this quantity in equation (3) is 
not unduly large. 

The effective mass coefficient k * can be obtained from 

IS* = k - Ar 


where 


/ o 

Ak = 2 / xK(x)dx 
^ 0 

Up to Xq ~ 0.4 the integrand plots approximately as a straight line 
against x, This fact allows an analytical solution 

Ak X Xo%(xo) (4) 

which is sufficiently acciorate for applications where Xq does not 
exceed 0.4. 
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CALCULATION OF PROPELLER CHARACTERISTICS 


For ideal conditions with hoth propellers ahsorhing equal power at 
equal rotational speeds and considering hoth propellers as being very 
close together, equations for determining the element loading coeffi- 
cients of dual -rotating propellers are given in reference 5 as follows: 


and 



i km 

1 + ^ KW sin'^0Q 



j (1 + w)w sin 
1 + I KW Bin^iZio 


K(x) 


( 5 ) 


(6) 


For the same ideal conditions, the aerodynamic helix angles are 
determined from 



where 




(9) 


As noted in the preceding discussion on the influence of a spinner on 
the K function, the k function in these equations should he replaced 
hy K*. Values of k and K(x) as functions of the parameter 

^ = J(l + w) are given in reference 5- These values for an eight - 

hlade dual -rotating propeller (including k ' for Xq = O.36) are given 
in figure 1. 


When the element loading coefficients and the helix angles have 
heen determined (equations (5) to ( 8 )), the element force coefficients 
can he calculated, provided the drag-lift ratio is known, from the fol- 
lowing equations: 
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= I 






cos ^ 
sln^^o 

cos ^ 
sln^0o 


(tan 


(tan 



^ + tsua 



( 10 ) 


(11) 


*^^r\ 2 (^ “ ^ ^f) 

dx /p “ X ^tan ^ + tan 7^ \dx 


' dC ^\ 2 (1 " ■*^“1 5 % ^ r )/^^ q \ 

^^Jb. ~ X ^tan + tan 7 ^^ ^^dx 


(13) 


The derivation of equations (lO) to (13) is given in the appendix. For 
the ideal conditions. 


« 


Wp = Vg = Constant 
K ' = Constant 

For nonideal cases, experience has shown that satisfactorily 
accurate calculation .of the characteristics of single -rotating propellers 
can he achieved hy considering the circulation functions to he dependent 
only on the local, value of w. The same assumption is applied for dual- 
rotating propellers. Dual rotation, however, involves another complica- 
tion, namely, the mutual Interference between the two propellers. In 
equations (7) and (8), for exaarple, the aerodynamic helix angles of the 
front and rear propellers axe different and this difference depends on 
the mass coefficient, which is unknown for nonideal cases. 

In the present calculations, for lack of better information, it is 
assiuaed that k' for the ideal and nonideal case is not greatly dif- 
ferent . At each station k ' is evaluated for the locsil value of w 
as though the whole propeller operates at the local value. 

Before preceding with the calculations, the propeller geometry, the 
fll^t Mach number, and the operating J must be given, and a set of 
suitable airfoil data must be available. When this information is given, 
the calculations are made by computing 

and tan (equations ( 5 ) to ( 8 )) at appropriate radial stations for 
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several values of v. The solidity a helng kaown., a value of cj is 
obtained for each v and the corresponding a Is determined from the 
airfoil data. The value of a for a given cj depends on the section 
Mach number which Is given by 


+ (Sf 

When oj' and ctj^ corresponding to each assumed value of w are obtained, 
compute 

and 

(Pr)^ = (5^ + 


These values are plotted against Vj the Intersection of these curves 
with the Pj. and given by the propeller geometry gives the 

operating values of and The operating values of Wp and W|^ 

are then used in equations (5) to (8) to compute the operating t 

and ^ for use In equations (10) to (13) (these quantities 

may also be determined by plotting the trial values against w and 
reading the values corresponding to the operating v-p and ^). Since 

the operating ccj and a of each propeller are now known, Cj can be 

determined and can be used to find the corresponding tan 7 from the 
airfoil data for the appropriate section Mach number. 


The power and thrust coefficients axe given by 



(14) 


dx 


(15) 
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The efficiency Is 


_ J(CTp + CTr) 

^ “ CPp + Cpj,- 


( 16 ) 


(IT) 


( 18 ) 


Reference has already "been made to figure 1, irtiich shows the K(x) 
and K values as functions of the parameter J(l + V) . The present 
calculations, in most cases. Involve unusually large values of J with 
the result tliat values of J(l + w) are heyond the range covered in 
reference 1 and extrapolation, as indicated In the figure, is necessary. 


FROraLLER AND AIRFOIL DATA 
Description of Propeller 


Blade-form curves for the NACA 3-(3)(05)-05 ei^t-blade dual- 
rotating propeller are given in figure- 2. NACA l 6 -serles airfoil sec- 
tiohs are used throu^out. The spinner -propeller radius ratio Xq for 
the particular configuration tested was O. 36 . Details of the test setup 
are given in reference 6 . The propeller was tested at forward Mach num- 
bers varying from about 0.35 to O. 925 . Blade-angle settings of the 
front propellers were varied from 55° to 80°. The rear blade angles 
were set at slightly smaller values, the object being to have each pro- 
peller absorb equal power at peak efficiency. 


The NACA 3- (3) (05) -05 propeller incorporates a pitch distribution 
which makes it suitable for operation at an unusually hl^ value of the 
advance ratio, about 7*2. The propeller consists essentially of two 


single -rotating propellers of opposite rotation. The 




curve 


of this propeller is shown for ccmparlson with the average curve for an 
ideal (Theodorsen) dual-rotating propeller with J = 7*15 iii figure 3* 
This comparison shows that, in terms of the loading at the O .7 radius. 
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the test propeller incorporates hl^er loadings outboard lower 
loadings inboard than does the ideal propeller of Theodor sen. 


Airfoil Data 

Airfoil data used in the calculations are shown in figures 4 to 9- 
Figures 4 to 6 show the variation of (or tan 7 ) as a function 

of section Mach number with cj as parameter. Figures 7 "to 9 shov the 
variation of a with M^ with cj as parameter. 

Data for these figures were obtained from reference 7* The test 
program of reference 7 involved the measurement of chordwise pressure 
distributions on a blade of the NACA 10-(3)(066) -033 propeller under 
operating conditions. By using the circulation functions for single- 
rotating propellers to determine the Interference velocities, it was 
possible to calculate the two-dimensional airfoil characteristics. 

Since the method of measurement gives the lift and the pressure drag 
only, it is necessary to estimate the friction drag. In the present 
case a friction drag coefficient of 0.004 was used. 

The propeller from 'which the airfoil data -were obtained has the 
same design lift coefficient as the propeller being investigated. The 
thickness distribution of "the two propellers differs, however, -which 
necessitates cross plotting against thickness ratio -fco obtain data for 
the proper thickness ratio. 

As previously mentioned, comparisons of airfoil data from the 
results of the special propeller tests -with standard two-dimensional 
wind-tunnel data show con^iaratl-vely small differences in the subcritical 
speed range. Beyond this range the differences are more pronoionced. ' 
Wind-tunnel data in the supercritical range are subject to largely 
unknown wind-tunnel--wall effects, and it is believed preferable -to use 
the propeller test data in propeller calculations. 


RESULTS AND DISCUSSION 


Calculations of the propeller characteristics -were made for 
several blade -angle settings and flight Mach numbers. The calculated 
total -power coefficients and efficiences are compared -with experimental 
values in figures 10 to 12. Figures I 3 to 15 give comparisons of the 
individual front and rear power coefficients. The blade-angle settings, 
the flight Mach numbers, and the values of J are indicated in the 
figures. 
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The calculated total -power coefficients are in very good agreement 
with the experimental values (figs. 10 to 12). This good agreement 
indicates that, on the average-^ the dueil -rotating propeller circulation 
functions give the interference velocities correctly and that the air- 
foil data show the relation between cj, a, and Mx properly. Com- 
parisons of the calculated individual front and rear power coefficients 
with the experimental values (figs. 13 to 15 ) show that the calculated 
values for the front propeller are always higher and the calculated 
values for the rear propeller are always lower than the experimental 
values. This discrepancy is believed to be due to an assumption in the 
derivation of the equations that both propellers operate with their 
planes of rotation very close together. 

Comparisons of calculated and experimental efficiency show good 
agreement between the two sets of data at M = 0.53 (fig. 10), but at 
the higher Mach numbers (M = 0.8 and M = O. 9 , figs. 11 and 12) the 
calculated efficiency is from 3 to 6 percent lower than the experimental 
efficiency. This discrepancy is probably caused by inaccuracies in the 
drag coeff iclenta" used in the calculations. Inaccuracies in the drag 
data would not appreciably affect the calculation of the power coeffi- 
cients because the drag vector in the torque direction is very small. 


CONCLUSIONS 


Theodorsen's propeller theory has been applied to the calculation 
of the characteristics of a dual -rotating propeller having nonideal 
load distribution. Calculations of the characteristics Involved the 
use of airfoil-section data obtained from the results of propeller tests. 
Within the range of the investigation, comparison of the calculated 
results with experimental results leads to the following conclusions. 

1. The calculated and experimental total -power coefficients are in 
good agreeiMnt at fll^t Mach numbers varying from 0.53 to 0 . 90 . 

2. The calculated and experimental efficiencies are in close agree- 
ment at a flight Mach number of 0.53- At flight Mach numbers of O. 8 O 
and 0 . 90 # the calculated efficiencies are from 3 to 6 percent lower than 
the experimental values. The discrepancy is believed to be due to 
Inaccuracies in the drag data. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APEEaroix 

DERIVATION OF FORMULAS FOR CALCULATION OF FORCE COEFFICIENTS 


The element torque Is 


<iQ 


With 


X 


r 

R 


and 


or 


Bh 

2itr 


tan 7 = — 


3 

^ pW^dCj COB 0(tan 0 + tan 7)x^ 

From reference 5 

and 

"k = 

'^O 

Inserting these equations for the resultant -velocities and di-viding 

hy pn^D^ gi-ves -yie Indl-vidual front and rear element torque 
coefficients 
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cos 0p 
sln^Q 


^an. ^ + tan 7p>jx^ 


(S^)r = I + f !^ + tan 7 r) x' 


The element thrust is given by 


dT = ^ pW^Bb(cj cos 0 - c^ sin 0^dr 


p II 

Substituting as before and dividing by pn D gives 


* I “ 'f"p * 


^)r = ■" 4 Vr - tan 0R tan 7 r) x 


It is readily seen that 


/ dCrj\ 

_ 2 

(t 

V^/F 

X 

(■ 


dx 


and 


'dCr 

dx 


2 {] - ’’b) /^a\ 

X ^tan 0 j^ + tan 7j^ J \^dx jj 
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(a) Mass coefficients. Xg for k' curve is O.36. 

Figure 1 .- Mass coefficients and circulation functions for eight-blade dual 
rotating propeller. (Short-dash line indicates extrapolation.) 





(b) Circulation functions 


Figure 1.- Concluded. 









Figure 2.- Blade-form curves for HACA 3“(3)(05)-03 dual-rotating propeller. 
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Figure 3»- Comparison of design blad.e loading of NACA 3-(3)(05)-05 dual- 
rotating propeller with "blade loading of optimum dual-rotating propeller 
J, 7.15. 



20 


KACA RM L51A17 


tan if' 



Figure 4.- Drag-lift ratio vsirlatioii used for station x = 0.45* 
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Figure 5-- Drag-lift ratio variation used for station, x = 0.60. 
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Saotlon Maoh number, 1^ 


Figure 6.- Drag-lift ratio variation used for stations x = 0.70 to 

X = 0.95- 
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Section Mach nuaber, 

Figure T*- Angle-of-attack variation used for station x = 0.45. 
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Figure 8.- Angle-of- attack variation uBed for. station x = 0.60. 
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Figure 9 -- Angle-of -attack variation used for stations x = O.7O to 

X = 0.95. 
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Figure 10.- Ccm^ariBon of eaperlmental t] and C] 

and Cp. M = 0.53* 


with ceilculated t) 






Figure 11.- Ccmpaxison of experimental ti and Cp with calculated ti 
and Cp. M = 0.80; Pp = 75°; Pr = 73°. 
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Figure 13 •- CCOTparlsou of experimeiital C„ and C„ with calculated 

Cp and Cp . M = 0 . 53 * 
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Adranoe ratio, J 

Figure lit.- Comparison of experimental Cp^ and Cp^ with calculated 
Cpp and Cp^. M = 0.80; % = 75°; 3p = 73°. 
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Figure 15 •- Conrparison of experimental C-p and Cp with, calculated 

•^F -^E 

Cp and Cp . M = 0.90j Pp = 75°; = 73°- 
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